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1. Navier-Stokes Equations: Background 

The present paper is an attempt to provide a brief orientation to Navier-Stokes 
equations from a probabilistic perspective developed in the course of working 
with a focussed research group in this area during the past few years at Oregon 
State University (OSU). An effort is made to furnish a bit more background for 
the uninitiated on some of the basics of these equations along with a summary 
description of the probabilistic methods and results. The approach is in some- 
what "broad strokes". The reader should be able to follow and/or supply most 
basic calculations, but the more technical proofs of some of the main results are 
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left to original references. More detailed accounts of some of the mathematical 
intricacies, conjectures and open problems identified by the OSU group can be 
found in Bhattacharya et al. [H||7j, Chen et al. \T2\, and m a recently completed 
PhD thesis by Orum 

A purely mathematical description of the Navier-Stokes equations naturally 
depends on the framework (function space) in which initial data and exter- 
nal forcings are provided, and in which solutions are sought, e.g. see see La- 
dyzhenskaya [6'2i \'6'6\ . Temam However, rather than begin with the tech- 
nical nomenclature, let us first take a look at the underlying physical picture. 
The 3-dimensional incompressible Navier-Stokes (NS) equations governing fluid 
velocities u(x, t) — (iii(x, i), u 2 (x, f), u 3 (x, £)), x e R 3 , t > 0, with initial data 
u(x, + ) = uo(x), are a manifestation of Newton's law of motion and mass 
conservation for a homogeneous (constant density) fluid: 



tions. The term p(x, t) is the (scalar) pressure, g(x, t) represents external forc- 
ing, and v > is the kinematic viscosity. Probabilists might want to write 



(4irvt) 2 exp{— -^L-} more familiar in applied mathematics. The incompressibil- 
ity condition refers to the requirement that 



The (in-viscid) equations corresponding to v = are referred to as the Euler 
equations, and are generally considered apart from Navier-Stokes. 

The equations comprise a (non-linear) system of four partial differen- 
tial equations governing four unknown scalar quantities; three coordinates of 
velocity ui, 1*2,1*3 and a scalar pressure p. The equations model the motion of 
a viscous incompressible homogeneous fluid in an inertial frame in free space. 
They are derived from conservation of mass and momenta of "fluid elements" 
(in the sense of continuum mechanics) under the additional assumption that 
there is a local linear relation between the stress and strain rates experienced 
by fluid parcels. In Q), the left-hand side + (u ■ V)u represents the accel- 
eration in an Eulerian frame of reference; i.e. one views the motion as that 
of a "parade of fluid elements" and u(x, t) is the velocity of the fluid par- 
cel at that location and at that time. Thus if a parcel has moved to location 
x(t) = (x(t),y(t), z(t)) at time t, then the velocity of the parcel is given by 
u(x,t) = (u 1 (x,y,z,t),u 2 (x,y,z,t),u 3 (x,y,z,t)) = (ff,^f,^f) and, using the 
chain rule, the acceleration is therefore given by the so-called material or con- 
vective derivative 




the term i/A as 



}^lvA in anticipation of the heat kernel in the form fc(x, t) 



V • u = 0. 



(2) 




(3) 
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The right side of £Q) resolves the forces acting on the fluid parcel. The first term 
yAu represents the viscous friction forces and is often referred to as the dissi- 
pative term. This term arises in the fluid mechanical derivation as the result of 
the aforementioned ad-hoc model of linearized shear stresses due to interactions 
between molecules comprising the fluid. The term g = g(x,y,z,t) represents 
body forces, such as gravity and/or other externally applied force fields. The 
pressure gradient term — Vp defines the local pressure force on a fluid parcel 
(depending on p up to an arbitrary constant). Physically the pressure gradient 
is a force which acts to maintain a homogeneous fluid density under incom- 
pressibility. Thus one may reasonably expect that pressure can be computed 
from the velocity and boundary conditions (at infinity in this case). We will 
see that the following orthogonality opens a way to mathematically first remove 
the pressure for focus on the incompressible velocity. Integrating the following 
"product formula" over a ball Bq (R) of radius R 

V • (pu) = u ■ Vp + (V • u)p (4) 

and then applying Gauss' divergence theorem yields the following version of an 
"integration by parts" type formula 

/ u ■ Vp = pun- / (V • u)p. (5) 

Jb {R) ' JdB (R) Jbo(R) 

Thus, if u is incompressible then the second integral on the right side of (JSJ van- 
ishes and the first integral is O(||pul(i?o(-R))||oo-R 2 ); here 1(Bq(R)) denotes the 
indicator function. Assuming pu decays sufficiently fast at infinity, one therefore 
obtains the following basic L 2 — orthogonality 

/ u-Vp = 0. (6) 

Remark. Decay assumptions are implicit to most of the theory to be described 
throughout this article, from orthogonality conditions ijBJ to Fourier transform 
formulations. Such boundary conditions at infinity on Navier-Stokes are the sub- 
ject of extensive investigation for both free-space and exterior domain problems 
since Finn's conception of the notion of "physically reasonable" solution; e.g. 
see Finn (T§|, Hey wood Galdi [231 - A perhaps more familiar illustration of 
the need to impose decay conditions occurs for example in the simple case of 
the heat equation, where one observes that since the Gaussian heat kernel is a 
tempered function, its convolution with a tempered distribution is well-defined 
as a function. In particular, it is sufficient to assume that the initial data is 
a tempered distribution for the well-known solution as a convolution with the 
heat kernel to apply. 

In summary, two distinguishing features of the equations JU are: (a) the non- 
linear term, referred to as the convective term, (u • V)u on the left-hand side, 
whose role in representing the acceleration is intrinsic and cannot be dropped 
by changing the model of the forces on the right-hand side, and (b) the incom- 
pressibility condition V • u = 0. This latter condition is an additional equation 
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to momentum balance obtained by consideration of mass conservation for a 
homogeneous fluid of constant fluid density. 

Systematic derivations of the Navier-Stokes equations can be found in any 
number of standard texts on fluids, e.g. Batchclor [3], Landau and Lifschitz 
[34] . Chorin and Marsden [15] . Thumbnail historical sketches of Claude Louis 
Marie Henri Navier (1785-1836) and George Gabriel Stokes (1819-1903) are 
available on the web. A brief textual history as it pertains to fluid flow as 
well as a fundamentally significant role in engineering design can be found in 
Anderson £Q. 

Let us briefly turn to the mathematical framework in which one might view 
these equations. Cannone |Hj contains an excellent survey of a variety of natural 
function spaces in which to analyze the equations. However, rather than attempt 
to develop the precise function space requirements at the outset for the results to 
be discussed, we will for the most part continue to confine ourselves to L 2 (R 3 ) 3 
and/or spaces of solenoidal, i.e. V • u = 0, tempered distributions u, where the 
reader will be able to check basic statements and calculations using standard 
Fourier methods as covered, for example, in Folland HJj|. The actual function 
spaces will evolve naturally in the course of the development of the probabilistic 
approach. 

An in-depth analysis of these equations given by Jean Leray [36] 1 has become 
a lasting benchmark for contemporary research. Under the assumptions that 
Uo,g € L 2 (R 3 ) 3 , V ■ Uo = 0, and an energy inequality of the form 

l/2||u(.,t)||l + f ||V • u(;s)\\jd S < l/2||u(.,0)|| 2 + f < g(-, S ),u(, S ) > ds 
Jo Jo 

(7) 

applies (in the distributional sense), Leray |3U proves the existence of a global 
(weak) solution to satisfying the inequality 10 . The corresponding unique- 
ness in L 2 (R 3 ) 3 of Leray 's solution remains open to this day for the case of n = 3 
space dimensions. From physical considerations, namely conservation of energy, 
the energy inequality 0) is expected to be an equality. However it has not been 
proven that such a solution exists in three-dimensions. A formal derivation of 
conservation of energy is obtained by taking the scalar product of the terms in 
with u and then integrating over R 3 . Assuming sufficiently rapid decay of 
u,p at infinity, integration by parts and incompressibility yield that the convec- 
tive term and pressure both vanish. For the scalar product with the dissipative 
term one has similarly by performing integration by parts 




= -H|V-u(,t)|| 2 . (8) 



Noting that J R3 |u(-,t)| 2 = J ^ • u one has the equation for conservation 
of energy 

~\H-,tm + v\\V ■ <;t)\\i =< g," > . (9) 

1 Robert Terrell has made his translation of Leray 's seminal paper to English available on 
the web at http://www.math.cornell.edu/~bterrell/leray.shtml 
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At another extreme for solution spaces, in defining the Millenium Prize Prob- 
lem for the Clay Institute, Fefferman ^7] stipulates that a physically reasonable 
solution must be infinitely smooth and decrease more rapidly than any power 
of the spatial variable x. However this choice of solution space is not with- 
out controversy in the international mathematical community. In a critique of 
this requirement, Ladyszhenskya remarks that "...one must leave both the 
choice of the phase space and the class of generalized solutions to the researcher 
without prescribing to him infinite smoothness or some other smoothness of 
solutions. The only requirement needed is indeed that the uniqueness theorem 
must hold in the chosen class of generalized solutions." Here she is addressing 
her basic point that the first fundamental question is whether the equations 
(and boundary conditions) do uniquely determine the dynamics of an incom- 
pressible fluid. In particular, questions regarding the qualitative behavior and 
properties of the solution naturally follow. In the context of this article bound- 
ary conditions refer to decay properties at infinity, however the methods extend 
to periodic boundary as well. 2 

Apart from definitions of classical and various forms of weak (in distributional 
sense) solutions, some of the mathematical intricacies are further illustrated with 
the work of Fujita and Kato Kato [301 m which a notion of "mild solution" 
is advanced as follows: The so-called Leray projection operator V is most easily 
defined by the orthogonal projection of the Hilbert space L 2 (R 3 ) 3 onto the 
closed subspace {u e L 2 (R 3 ) 3 : V ■ u = 0}. In particular, denoting (component- 
wise) Fourier transforms on u = (iti,ti2,«3) G L 2 (R 3 ) 3 by u = (ui, 112, U3), 
where for / € L 2 (R 3 ) 



the subspace of divergence-free (incompressible) vector fields corresponds to the 
following basic orthogonality in Fourier frequency space 



Now, since orthogonality of the projection requires that u(£_) — Vu(^) be orthog- 
onal to T'u(^), this vector difference is in the direction of £, i.e. u(£) — Vu^) — 



(if, ■ u(0)|||- Thus, 

^(0 = fi(0-(ilr-u(0),l| =(i-i^)u(0. £eR 3 , (12) 



where I = (($ij))i<i,j<3 is the identity matrix and the tensor product (£> of 
two vectors a, b £ R 3 is defined by the matrix a ® b = ((cii&j))i<i,i<3- The 

2 For the most part analysts do not make a significant distinction in the choice between 
free-space and periodic boundary conditions. While both free-space and periodic boundary 
are amenable to stochastic cascade representations, there are subtle differences due to the role 
of small non-zero frequencies which do not arise in the wave-number lattice associated with 
Fourier series. 




(10) 



£ • Pu(0 = 0. 



(11) 
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projection extends to other function spaces for which Fourier transforms can be 
defined. In view of the incompressibility condition and one may apply the 
projection operator to Q to obtain for incompressible initial data Uo belonging 
to a suitable function space, e.g. L 2 (R 3 ) 3 , the "projected equation" 

— +P(u-Vu) = ^Au + P(g), u(x,0+) = u o (x). (13) 

If one views 11130 as a (vector) heat equation forced by the (unknown) term 
w = V(g — u • Vu) then a standard Duhamel argument, e.g. see Bhattacharya 
and Waymire 4 , may be applied to deduce the following integral equation 

u(x,i) =T t u (x)+ f T t _ s w(x,s)ds, (14) 
Jo 

where {Tt = e tl/A : t > 0} is the (self-adjoint) semi-group associated with 3- 
dimcnsional Brownian motion Z = {Z(t) : t > 0} with zero drift and diffusion 
coefficient 2v, i.e. T t g(x) = E x g(Z(t)) = / R3 g(y) 1 4 cxp{- 3 ^ ? |y-x| 2 }dy, 

and applied component- wise. Solutions to 1)14(1 satisfying V ■ u = are referred 
to as mild solutions. This is the starting point of Kato's approach in which he 
proved existence of mild solutions for "small" initial data in the special function 
space L 3 (R 3 ) 3 . Specifically, ignoring external forcing for simplicity, Kato [3Hj 
showed that there is a number A > such that for ||uo||l 3 < A one nas a 
solution such that t — * u(-,t) G L 3 (R 3 ) 3 is continuous and satisfies ((14(1 on 
every compact time interval [0, T], T > 0. However it wasn't until very recently 
that uniqueness was proven by Furioli, Lemarie-Rieusett and Terraneo |22| . 

The Picard iteration scheme provides a useful technique naturally associated 
with the (projected) Navier-Stokes equation. Observe that using incompress- 
ibility the non- linear term may be expressed, analogously to uu x — ^(u 2 ) x in 
one-dimension, in the equivalent form 

(u-V)u = V-(u®u). (15) 

Then the iteration is given by 

u n+ i(x,t) = F(x,t) - &(u„,u„)(x,£) (16) 

where F(x, t) =T t u (x) + /„* T s Pg(x, t- s)ds, 6(u,v) := /* VT t - s V ■ (u® v)ds, 
u(°l(x,t) =T t u (x) and u x (x,t) ^ F(x,t) + b(u^ ,u^){x,t). The convergence 
of the iterates follows from showing that that these are iterates of a contraction 
on a ball in a suitable function space; c.f. Meyer (2004) and references therein. 

It may be noted that the function space L 3 (R 3 ), or more generally L"(R"), 
enjoys the very special norm scale invariance of the form 

IIMU» = /A(x) = i/(Ax),xeR n ,A>0. (17) 

Such Banach spaces are referred to as limit spaces. Excellent surveys of contem- 
porary literature which grew out of Kato's approach are given in Meyer gU, 
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Cannone ^3 EJ ■ A more general notion of limit space arises in the context of 
majorizing kernels and will be given in Section|3] see (|64f> . The Picard iteration 
approach may also be viewed as a starting point for the probabilistic methods 
initiated by LeJan and Sznitman |37j . though perhaps not in a way one might 
expect on first glance. Their's is based on a branching representation in Fourier 
space in which the non-linearity in the Fourier-transformed version of is as- 
sociated with a "binary branching product" of certain random walks in Fourier 
frequency space. The essential difference in the two approaches is in the nature 
of bounds sought for convergence of two distinct but closely related limits. In 
the former one seeks contraction bounds for convergence of the iteration scheme, 
while in the latter one seeks integrability bounds on stochastic variables, defined 
by a.s. limits of a stochastic iteration, to ensure finiteness of their expected val- 
ues. To highlight the close connections between these two approaches it may be 
noted that the aforementioned uniqueness proof by Furioli, Lemarie-Rieusett 
and Terraneo [52] for Kato's solution in L 3 (R 3 ) was obtained by arguments, 
specifically estimates on the associated bilinear form b(u, v), which the authors 
explicitly acknowledge to have been inspired by the probabilistic approach to 
uniqueness for a certain other space obtained by LeJan and Sznitman |37|: 
namely the space T = {u e <S'(R 3 ) : sup ? |£| 2 |u(£)l < °°} or > more accurately, 
the solenoidal distributions belonging to this space T . A distribution whose 
Fourier transform is essentially bounded is referred to as a pseudo-measure in 
harmonic analysis. Thus an alternative description of this space is one of locally 
integrable tempered distributions for which Au is a pseudo- measure; now often 
denoted VM 2 where the superscript 2 refers to the |£| 2 . pre-factor, c.f. Cannone 
JS]. In fact this may also be shown to be a Besov space; see Cannone and Plan- 
chon For the unitiated, a brief historical orientation to the various such 
function spaces which arise in this context may be found in Peetre |48| . 

Since, loosely speaking, algebraic products transform to convolution products 
of Fourier transforms and since derivatives transform by simple algebraic factors, 
the transformation of the non-linearity under Fourier transform in the (quasi- 
linear) Navier-Stokes equation is deceptively similar to the convolution term 
one obtains in Fourier space for (semi-linear ) reaction-diffusion equations, e.g. 
KPP. But one must not forget about the incompressibility condition, which is 
completely absent for the case of reaction-diffusion. Nonetheless, it is natural to 
search for branching Brownian motions in physical space; e.g. see McKean jJU] 
in the case of KPP. In the case of the scalar Burgers equation, again devoid of a 
notion of incompressibility, one may identify such a representation; see Waymire 
[5fi| . Ramirez |49j . Moreover, recently a probabilistic representation involving 
scmi-Markovian branching Brownian motions in physical space was obtained by 
Ossiander 07] in which some basic notions from probabilistic potential theory 
make a natural appearance. In the next section we will further indicate natural 
probabilistic ties to harmonic analysis in the way of the "background radiation" 
process in physical space introduced in Gundy and Varapoulos j27J and Gundy 
and Silverstein |2f>| . 

The previously noted basic results capture the essence of the paradigm fur- 
nished by iterative methods in a wide range of function spaces. Namely, one 
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either obtains for sufficiently small initial data Uo existence and uniqueness of 
global solutions, i.e. for all times T > 0, or one proves existence and unique- 
ness for arbitrarily large initial data Uo for short durations T = T(uo) < oo. 
A concisely formulated point for caution concerning semi-group approaches was 
provided in a paper by Montgomery-Smith 03] which suggests a much deeper 
role for incompressibility than has thus far been captured by these methods. 
However there are many significant questions beyond existence and uniqueness, 
e.g. numerical computation and Monte-Carlo simulation algorithms, boundary 
value problems, and steady-states, which stand to benefit by continued devel- 
opment of probabilistic methods. It is manifestly clear that the introduction 
of probabilistic reasoning has certainly lead to new insights into the structure 
of these equations and their solutions. On the other hand, it is also clear that 
incompressibility is being undervalued in these analyses to date. 

2. Vorticity, Pressure, Incompressibility and Background Radiation 

One need do no more than point to the classic papers on diffusions by W. 
Feller, to name just one example, to highlight and exemplify the impact that 
the analytic theory of linear (parabolic) partial differential equations has had on 
the development of probability theory. On the other hand there is no question 
to the significant counterpoint that the subsequent development of martingale 
theory, the Ito calculus and the Malliavin calculus has provided in the study of 
such pde's. However, unlike the fundamental equations governing diffusion, the 
historic origins of the Navier-Stokes equations are devoid of explicit probabilistic 
reasoning for their formulation - making the appearance of probabilistic methods 
here all the more mysterious from the outset. 

There is that Laplacian term! One of the earliest applications of probabilis- 
tic ideas in connection with incompressible Navier-Stokes was introduced by 
Chorin |13j . based on the associated vorticity equation. The vector cross prod- 
uct w(x, t) = V A u(x, t), i.e. the curl of the velocity field, defines the vorticity. 
Although the physical equations are derived in n — 3 dimensions the formula- 
tion extends to n > 2 dimensions. In n . — 1 dimension the notion of vorticity 
is lost; the corresponding scalar equation for velocity ^ = + + g 
being referred to as Burgers equation. However in the lowest dimensional case 
(n = 2) of velocity u = (ui, U2) in fll}. one may also identify the vector vorticity 
u> = (0, 0, lu) with the scalar component u). 

Gradients are annihilated by taking the curl since for any sufficiently smooth 
scalar function a 

VAVfl = V A (da/dxi, da/dx2, da/dx 3 ) 

— (d 2 a/dx 2 dx 3 — d 2 a/dx 3 dx 2 ,d 2 a/dx 3 dxi — d 2 a/dx\dx 3 , 

d 2 a/dx\dx2 — d 2 a/dx2dxi). 
= 0. (18) 

In particular taking the curl provides another approach to remove the (unknown) 
pressure gradient in QJ. Applying some simple vector identities, e.g. see Chorin 
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and Marsden [El pp. 163-164], namely (u ■ V)u = ±V(|u| 2 ) + uj A u and V A 
(wAu) = w(V • u) — u(V • ui) + (u • V)cJ — (uj ■ V)u, and using incompressibility 
and (|18fl . one may express the curl of u satisfying (JTJ as 

^ = ^Acj - (u • V)co + (uj ■ V)u + V A g, V-w = 0. (19) 

A solution uj of l|19|) may then be used to obtain a solution u to Q and, in turn, 
the pressure p as follows: According to the Helmholtz- Hodge decomposition any 
smooth vector field u which decays sufficiently fast at infinity may be uniquely 
represented as the superposition of a gradient and a (necessarily incompressible) 
curl 

u = Ul +u 2 , (20) 

where 

Ui = -V$, u 2 = VA* (21) 

for a scalar potential $ and vector potential W obtained by solving a (vector) 
Poisson equation 

A$ = -Vu, A* = -AAu, V-# = 0. (22) 

In particular one may recover u from vorticity uj via the following so-called 
Biot-Savart law. 

u(x,t) = L*w(pc,t), (23) 

where L is a linear (convolution) kernel furnished by the divergence free solution 
^ to the Poisson equation, i.e. by inverting the Laplacian via (for n = 3) 

*(x,t)-A- 1 (-c) = --L f ^1 (;M ) dVj v* = 0, (24) 

47r Jw> l x - y| 

and hence (for n = 3) 

u(x,i) - VA (-1^/ y^—Mt,y)dy)=^-{ V x -^—oj(t,y)dy 
4tt 7 r3 |x - y| 4tt 7 r3 |x - y| 

= £*w(x,t). (25) 

Similarly for one needs only to apply the appropriate potential kernel to 

invert the Laplacian. In particular for n — 2 dimensions one has 

u(x, t) = VA(-i / logflx - y|)w(y, t)dy) = L * oj{t, x), (26) 

^ JR 2 

where 

M^rns^rii), x = (ajijara) € R 2 . (27) 

Z7T|X| Z 

For numerical computations Chorin |13j viewed (|19|) by an operator splitting, 
i.e. along the lines of a Trotter-Kato semi-group product, into successive incre- 
mental advective and diffusive flow processes using distinct numerical schemes. 
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In this connection, assuming n = 2 observe that the term w ■ V in (|19[l van- 
ishes for n = 2; here one should recall the convention of using the scalar oj 
and vector (0,0, w) interchangeably in n = 2 dimensions. Thus, assuming also 
g = 0, the two-dimensional vorticity equation describes the evolution of a scalar 
oj = w(x, t) by 

= i/Aw - (u • V)w. (28) 

Now, if u were given (computed) over an increment of time, then the scalar 
component of vorticity would appear to be governed by a diffusion equation. 
In this regard, denote two-dimensional standard Brownian motion by W = 
{(Wi(t),W 2 (t)) :t>0} and let Z := {(Z x (t), Z 2 {t)) : t > 0} be the diffusion 
defined by the associated stochastic differential equation 

dZ(t) = V2vdW(t) - u(Z(i), t)dt (29) 

Then one is naturally led to consider 

ui(x,t) := Eui (Z(t)). (30) 

This is the starting point for an approach which has been exploited heavily 
as a numerical Monte-Carlo method, and has been studied closely by analysts 
and probabilists alike; e.g. Chorin |T3], Marchioro and Pulverenti Goodman 
|2"1) . Long |3H], Szumbarski and Wald |FT], Meleard [H]. In this regard, Bus- 
nello, Flandoli, and Romito |S] also obtain a natural probabilistic version of the 
Biot-Savart law in three-dimensions. While we will not pursue these particular 
approaches any further, we will next take the idea of considering linearization 
and Duhamcls principle in an apparently new direction which exploits incom- 
pressibility in an interesting way. 

To recover the pressure from its projection onto a divergence free vector 
field let us recall i|12|) ■ The matrix elements ^ may be be viewed as Fourier 

symbols of iterated Riesz transforms RiRj, i.e. Rjf{£) = — note the 

signs and constants depend on the version (|I0|1 of Fourier transform used. In 
particular, on appropriate function spaces, e.g. L 2 (R 3 ) 3 , one has 

r = i + n, (31) 

where 1Z = ((RiRj))i<i.j<3 is the matrix of iterated Riesz transforms. To exploit 
incompressibility in JI} note that the terms linear in u are divergence free. Thus 
take the divergence to annihilate the linear terms d ^ t u = = A(V • u) = 
V • (Au), and assume no forcing g = 0, so that one has 

V • (u • V)u = -V • Vp = -Ap. (32) 

Taking Fourier transforms with continued use of incompressibility one arrives 
at the following determination of the pressure 



P= ^ RjRk{ujU k ). 

l<j",fe<3 



(33) 
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Apart from what may be gleaned from its potential theoretic role in making 
projections onto divergence free vector fields, a basic understanding of incom- 
pressibility remains illusive from a probabilistic perspective. However, even in 
the case of linearized Navier-Stokes equations, the incompressibility condition 
plays a fascinating role. Two common linearizations are (i.) Stokes equations and 
(ii.) Oseen equations, respectively given by linearization around zero velocity or 
by linearization around a non-zero constant U, respectively. Specifically one has, 
respectively, for Stokes 

" U ;Au-Vp + g, Vu = 0. (34) 



and for Oseen 

~dt 



dt 



U • Vu = i^Au - Vp + g, V ■ u = 0. (35) 



For simplicity and without any significant loss of generality let us consider the 
Stokes equation H34fl . Observe that for given twice continuously differentiable 
and incompressible initial data Uo = {uo~\ Uq 3 ^), one has that in the un- 
forced case (g = 0) the solution (u,p) = ((iti, u^, u^),p) is given by the familiar 
semi-group representation 

Uj (x,t) = EvtfW^Wit)), j = 1,2,3, p = c, (36) 

where W = {(Wj.(i), W<z(t), W^(t)) : t > 0} is three-dimensional standard Brow- 
nian motion, and c is an arbitrary constant. In particular the linearity in u and 
the apparent affine form of (|34J) due to — Vp are reconciled by the incompress- 
ibility. In other words, the incompressibility of u (at times t > 0) is directly 
inherited from that of Uo, causing — Vp to instantaneously "drop out" at t = 0. 
Accordingly one may view the matrix-valued Gaussian (heat) kernel 

K(z, t) = diag[A:(z, t),k(z, t),k(z, t)}, 

3 M 2 

k(z,t) = (47wi)^e~^ (37) 

furnished by the transition probabilities of the Brownian motion, as the funda- 
mental solution to l|34|) for incompressible initial data. 

Let us consider, on the other hand, the evolution of smooth but compress- 
ible initial data u . Also let us include a possibly compressible forcing term g. 
We seek an incompressible vector field u(x, t) = 7 5 u(x, t),t > 0, such that by 
orthogonality of Vp to incompressible u for t > 0, 

dVu 

— = vAVu + V(g), t>0, Pu(x,0+)=?u o (x). (38) 

The compressible initial data Uo will necessitate that u(x, t) = 7 5 u(x, t) for 
t > 0, be discontinuous at t = 0. Applying Duhamel's principle to (|38|l one 
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arrives at 

u(x,t) = / K(z - x,f)-Pu (z)dz + f [ K(x- z,t- s)Vg(z, s)dsdz 
= / r(z - x,t)u (z)dz + / / r(x- z,t- s)g(z, s)dsdz, (39) 

where 

r = {I + ft}(K), (40) 

and 7?. = ((RjRk))i<j,k<3 is the matrix of iterated Riesz transforms. The sim- 
plest way to check the commutativity required for the second equation in the 
display ((39(1 is by Fourier transforms along the lines indicated earlier for (|31|) . 
To compute the pressure one simply takes the divergence in (|34() and use in- 
compressibility to get 

V • Vp(x, t) = Ap = V • g, t > 0. (41) 

For £ = one may apply the Hclmholtz- Hodge decomposition to recover uq from 
■pu defining the pressure at t — by 

p(x,0) = $ = A" x (-V ■ u ). (42) 

In particular, therefore, 

p(x, t) = A" x (-V ■ uo)(x)5 (*) + A" X (V ■ g)(x, i) (43) 

One may note that 1)39(1 actually embodies two Duhamel principles; one applied 
and one derived! Namely, we apply a Duhamel principle for Brownian motion 
to (|38|l with "P(g) as the forcing, and derive a Duhamel principle for 134fl with 
arbitrary initial data using the fundamental solution (|40() of the unforced linear 
incompressible pressurized equation and adding g as the forcing. 

The representation 1(40(1 of the (signed) fundamental solution T of Stokes 
equation for (u,p) is the well-known classic formula of Oseen a l so see 

Solonnikov 50 , Thomann and Guenther |53| . One may check a semi-group 
property of the form 

r(x,t+s)=/ r(x- y ,t)r( y ,,s)dy. (44) 

Moreover, V has the following equivalent representation 

T = K + Hess(V>), Aip = -K, (45) 

where Hess denotes the Hessian matrix. A simple way to check this is via Fourier 
transforms since %fr is the Fourier symbol for the operator — A -1 . Thomann 
and Guenther |53j recently made an explicit computation of the full fundamental 
solution in terms of special functions which has had major implications for 



E.C. Waymire/ Probability & incompressible Navier-Stokes equations 



13 



analysis of Stokes equation in unrestricted domains. Their precise formulae may 
be expressed as 

A- 1 JT=- rT l-W a " xl 



4tt|z-x| y/tt 

15 |z-x' 2 
3 lFl(1 '2'^ 

(Xj - Zj)(x k - Zk) ( H 5 |z - X 



-RjRk(K) = K(z - x, *){-- ^(1, -, L—L)s jk 



^(1,-,^)-!^}, (46) 

where Erf(s) = fo e ~ y dy denotes the error function and iFi(a,c,u) — 

^2n=o (c)^ ^ 0T = r7E) > * s a connuen t hypergeometric function. This 
is a hugely significant formula in view of the rich mathematical literature on 
parametric dependencies, asymptotic growth and/or decay properties etc of the 
special functions which appear in the formula (|46|l . 

The matrix T(x, 0) is the orthogonal projection of the (vector) Dirac distri- 
bution (5 X onto divergence-free vector fields. Thus, the effect of incompressibility 
is instantaneous in that the solution is obtained by projecting the initial data 
onto divergence-free vector fields and then solving (|34f> with incompressible ini- 
tial data as above. In fact one may observe that 

r(x,i)=/ K(x-y,t)r(y,0)dy, x. — (x, y, z) G R 3 . (47) 

The transformation r(xo,0)uo is the (instantaneous) projection of Uo onto its 
divergence free (incompressible) component. 

The sign changes and discontinuity at t = in the fundamental solution 
make the issue perplexing from a probabilistic standpoint. However, there is 
some hope that a probabilistic unraveling of this semi-group might lead to ways 
to deal with boundaries. Thomann and Guenthcr (personal communication) 
have pointed out that if one poses a Neumann boundary condition on the half- 
space, the usual "reflection method" applied to the fundamental solution T for 
the unrestricted problem fails. 

In view of the prominent role of Riesz transforms it seems natural to explore 
a role for Gundy- Varapoulos-Silverstein's background radiation process in the 
present framework. Thus we close this section with elements of such a new 
connection. The background radiation process Z = {Z t = (Yf,X t ) : — oo < 
t < 0} belongs to the class of "approximate Markov processes" defined, as a 
collection of measurable functions, on a sigma-finite measure spaces (fi, J 7 , P) 
taking values in the half-space H3 = [0, 00) x R 3 . Foregoing its construction, 
the process is defined by the following properties (Gundy |25p : 

1 {Y t : —00 < t < 0} is distributed as a one-dimensional Brownian motion 
on (0,oo) started at (arbitrary) yo > and killed upon entry to the state 
0. 

2 lim t jo Yt = a.s. 
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3 {X t : — oo < t < 0} is an independent Brownian motion with initial 
distribution (at time H = — oo") given by Lebesgue measure on R 3 in the 
sense that for each y > 0, defining r y = inf{< < : Y t = y}, the process 
{(y, X Ty+t ) : < t < —T y } is distributed as a Brownian motion started 
with Lebesgue measure on the hyper-plane {y} x R 3 at t = 0. 

4 X = lim t | X t exists a.s. and is distributed as Lebesgue measure on R 3 , 
i.e. the induced measure P o Xq" 1 is Lebsesgue measure. 

These properties were used by Gundy and Varapoulos j2H in Gundy and 
Silverstein to prove the following formula for iterated Riesz transforms of a 
function / belonging to Schwartz space S 

RjR k f{x) = bE( [ A jk Vf ■ dZ t \Z a = (0,x)) (48) 

J — oo 

where /o denotes the harmonic extension of the function / on the R 3 to the 
half-space H3, the matrix Ajk = Bj ® e^, i.e. having 1 in location (j,k) and 
0's elsewhere, defines the indicated martingale transform. With a little more 
work one may in fact calculate b = — \ using their indicated methods. In view 
of l|4U|l and viewing conditional expectations as disintegration formulae, one 
may express the solution to Stokes linearization (|34l) in terms of background 
radiation as follows, writing u = (it^ 1 -*, u^ 2 \ u^): 



tt W) (x,t) 



/ fc(z-x,t)u[, j) (z)dz 

J-R? 

3 

+ ^2R ] R k (K){z-x,t)u {k) {z)dz 

k=l 

= E Xo=x u^(X* t ) 
3 

+ Y, ^(^' fe) (x,t)|Zo-(0,z))4 fe) (z)dz, (49) 

where {X£ : t > 0} is the reversed-time process (Brownian motion) introduced 
in Gundy and Silverstein |26j . and 6^ k ' is the martingale transform defined by 

e W = ^ fe )( x ,i) =-\ I ej ® e k Vif (Z s ) • dZ s (50) 
^ J-00 

and for fixed t, x, Kq is the harmonic extension of the the heat kernel z — * K(z — 
x, t) to the half-space H3. In particular, defining the matrix- valued function 

e(x,i) = ((^ fe )(x,t))) 1 < J , fc < 3 , (51) 

one has the representation 



u(x, t) = ^ Zo =(o,x)Uo(X t *) + £(6(x, *)uo(Z )). 



(52) 
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Form here one may regard the non-linear term in P as a forcing and exploit 
a Duhamel principle to obtain an integral equation involving background radi- 
ation. While this connection has not been pursued beyond this survey article, 
an alternative real-space semi-Markov cascade representation of Navier-Stokes 
is fully developed in Ossiander j^Zj which also contains intriguing connections 
with the fundamental solution obtained by Thomann and Guenther |53| for 
linearized flow. 

3. Multiplicative Stochastic Cascade 

One may give purely fluid mechanical derivations of diffusion equations based 
on mass conservation and linear flux laws such as Ficks and Fourier laws, how- 
ever the most compelling physical arguments are no doubt Einstein's implicit 
use of the central limit theorem and law of large numbers. While the modern 
approach of hydrodynamic scaling limits of interacting particle systems seeks to 
arrive at the Navier-Stokes equations from a more intrinsic probabilistic or sta- 
tistical mechanical foundation , such approaches are fundamentally much more 
complex than the connections that are possible between Brownian motion and 
the heat equation; see e.g. Esposito, Marra and Yau JH], Yau That said, 
owing to a remarkable insight of LeJan and Sznitman |H7j it is nevertheless pos- 
sible to approach certain classes of solutions to these equations by probabilistic 
constructions of branching random walks and multiplicative cascades in a way 
which preserves some basic Markovian semi-group structure. This is the topic 
to be briefly described in this section. 

As remarked earlier the multiplicative cascade approach appears as a math- 
ematical device, closely related to Kato's contraction map techniques for ob- 
taining a mild solution on appropriate function spaces. While a pure analyst 
may take solace in this realization, as a probabilistic technique it has additional 
virtues. Chief among these is the "natural" emergence of solution spaces (and 
their attendant norms) for which one may obtain unique global solutions under 
a "small" initial data constraint as a consequence of the probabilistic represen- 
tation. The calculation of steady state solutions are also quick to follow in these 
function spaces; see Bhattacharya et al. 0. New questions and ideas with regard 
to Monte Carlo methods emerge naturally in more pragmatic considerations of 
numerical solutions; see Ramirez |49|. However, whether this approach can lead 
to truly significant new results on the existence and uniqueness of global solu- 
tions under smooth initial data without regard to "size" remains to be seen. It 
is clear that a more substantial use of incompressibility is lacking in the existing 
theory; c.f. Montgomery-Smith |44|. 

Example (Warm-up). Some of the basic ideas under multiplicative cascade 
representations are more generally applicable to diverse classes of evolution 
equations, including certain linear parabolic and fractional diffusion equations, 
semi-linear reaction-diffusions and some quasi-linear equations such as two- 
dimensional incompressible Navier-Stokes equations and one-dimensional Burg- 
ers' equation; see Bhattacharya et al. Although interesting from a number 
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of other points of view, e.g. physical space representations, Monte-Carlo sim- 
ulations, etc., these examples are not capable of furthering understanding of 
incompressibility. That said,, the following extremely simple example serves to 
illustrate some of the most basic graph theoretic and probabilistic ideas involved 
in this approach. Consider 

u t — aAu + b ■ Vu u{x, 0) = uq(x), (53) 

in n > 1 dimensions, where a > 0, and b € R™ are constants. To quickly get the 
flavor of the method, consider the spatial Fourier transform of (|53[1 

fi(£, t) = u (0e- al(l2t + ^| /* a\ee- a ^ 2s u(Z, t - s)ds. (54) 

a Kr Jo 

Now consider the random linear tree Tg(t) rooted at a vertex 9 of type t;g — £ 
which, after an exponential length of time is replaced by a single vertex (1) of 
the same type = £. Proceeding in this manner one may calculate that the 
very familiar solution exp(— a|£| 2 + ib ■ £)uq(£) is the expectation of the random 
product X (6, t) initialized by £g — £ and consisting of factors m(£) = ib ■ £/a|£| 2 
at each vertex until termination where one attaches the end factor Uo(£)- The 
indicated stochastic recursion then yields 

X(0,t) = m(O" (t) fio(0, 
and the solution is represented as 

fifr t) = E (e=i X (9, t) - E is ^m(0 m M0, (55) 

where N(t) is the Poisson process with parameter A(£) :— a\£\ 2 which counts 
the number of times the exponential clocks ring before time t. In particular the 
Poisson process formally suggests a Fourier dual role to that played by the stan- 
dard Brownian motion in the real space expectation formula. Similarly one may 
obtain a formal "dual" Feynman-Kac formula (in the sense of this Fourier trans- 
form representation) under the complex measure condition on coefficients given 
by ltd [35] ; see Kolokoltsov [31] , and Chen et al. JHI ■ I n particular this approach 
makes Ito's complex measure condition completely natural from a probabilistic 
point of view. One may also obtain a Fourier version of McKean's |4()j branching 
Brownian motion formula for KPP, and other interesting equations; see Orum 
|45| . Also included are, for example, the generalized fractional Burgers equation 
of the type considered by Woyczynski, Biler, and Funaki 54 . 

The essential idea introduced by LeJan and Sznitman |37| and presented 
here in the generalized form developed in Bhattacharya et al. [H[Zj is as follows. 
Without loss of generality assume that the forcing g is incompressible; else 
replace g by Pg below. Either by directly taking Fourier transforms in the 
projected equation and then introducing an exponential integrating factor, 
or by taking Fourier transforms in l|14|l. one may deduce the following integral 
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equation in Fourier frequency space: 

u(£,i) = e-^l 2t u„(£) + /V"l^{|£|(27r)-i 
Jo 

/ 11(77, t - s) ® s u(£ -r),t- s)dr) + g(£,t - s)\ds, 
where for complex vectors w, z 



(FNS) 



to (gi^ z — — i(e£ ■ z)tt£±w, = j^r, and 7r^j_ui = w — (ej • w)e^ (55) 

is the projection of w onto the plane orthogonal to £, and ^ > is the viscosity 
parameter. For £ 7^ 0, LeJan and Sznitman [37] rescale the equation (|FNS|I to 
normalize the integrating factor e~^' 2s to the exponential probability density 
H£| 2e_ly '^' s ■ They then observe that the resulting equation is precisely the form 
for a branching random walk recursion for %(£, t) := i^|£| 2 u(£, t), for which the 
kernel |£— 7/|~ 2 |?7|~ 2 is naturally constrained by integrability to dimensions n > 3 
for normalization to a transition probability. Alternatively, Bhattacharya et al. 
introduce a Fourier multiplier 1/h, where h(£),£ € M^/i := £ R 3 : ^ 0} 
is a positive function such that < h * h(£) < 00, which is used to re-scale the 
Fourier transformed equation (|FNS|) by factors l/h(£) in compensation for the 
temporal normalization of the exponential factor to a probability. Namely, we 
consider the equation JFNSfcH obtained from (|FNS|) as 

X (U) = e-»W 2 XQ (0 + [ We-W'^m® [ xfoi,f-«) 



W h xW h 



x(m,t-s)H^,d Vl xdri 2 ) + ^,t-s)}ds, £&W h (FNSh) 



where 



2h*h(£) MO 2g(g,t) 

m(0 = K^M0' *° (0 = W ^' t) = HW (56) 

and H (£, c??7i x d^) is for £ € Wh the transition probability kernel, with support 
contained in the set { (771 , 772 ) € Wh x Wh : r]i + r/ 2 = £}, defined by 

/ f(r ) i,V2)H(Z,d m xdr l2 )= [ m ~ V,v) H [~ d V (57) 

for bounded, Borel measurable / : VF^ x Wh — ► R. In view of the conservation of 
wave-numbers 771 + 772 = £ reflected in the support of the transition probability 
kernel, the definition of the supporting set of frequencies Wh may be relaxed 
to a semi-group Wh Q {£ € R 3 : £ ^ 0}. However in general we include the 
following additional exterior condition in defining jFNS/J : 

*(£,*) = 0, £eW£,t>0. (58) 
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One may also note that with the exception of du/dt and g, the Fourier transform 
of each term in vanishes at zero frequency. Thus the value of u(0, t) = 
J R3 u(x, t)dx is explicitly determined from the initial data and forcing by 

^u(0,*)=g(0,t), u(0,0)=u (0). (59) 

The probabilistic approach is based upon an interpretation of the re-scaled in- 
tegral equation. This is achieved in terms of expectation values of multiplicative 
cascade solutions to stochastic recursions generated by a multi-type branching 
random walks in Fourier space. The process may be initiated at the root vertex 
9 of a binary tree graph having wave- number type £ ^ 0. After an exponential 
holding time Sg with rate v\£,\ 2 an independent coin flip is made. If the outcome 
is tail, denoted [k$ = 0] , then the tree terminates as a single edge of length Sg 
connecting 9 and < 1 >. If the outcome is head, denoted [Kg = 1], then the 
edge of length Sg branches into two vertices < 1 > , < 2 > of respective types 
£i,£2,£i + & = £, selected according to the probability kernel if (£,•). Given 
£ two independent exponential clocks with respective rates f|£i| 2 and v\&\ 2 
are reset at the vertices < 1 > , < 2 > and the process is repeated with each 
of these vertices as "roots" . The root 9 has an associated birth time Bg = 
and an end-time Sg at which it is terminated or it branches. Similarly each de- 
scendent vertex v =< v\, . . . , v k >,Vj £ {1, 2}, has a corresponding birth time 

B v := Y^ k j=i S<vx,...,Vj>- and an end-time B v + S v . At a given time t > one 
assigns multiplicative weights to vertices v such that either [B v + S v < t. k v — 0] 
or [B v < t < B v + S v }. In the first case the re-scaled forcing f(^ v ) is assigned, 
while in the latter case the re-scaled initial data x(£i>) is assigned, where 
denotes the frequency type of v; see Figure ^ 



One may now use the stochastic branching model represented by the col- 
lection of random variables {( v ,k v ,S v : v e V := U^ {l,2} fe } to recursively 
define a random functional related to (|FNS|I through its expected value. Namely, 
for measurable functions xo : Wh ~ > C 3 and <p : Wh x [0, oo) — > C 3 , and for 
£e = £ G Wh, t>0, the stochastic functional X (9, t) is recursively defined start- 
ing from v — 9 by 

)Xo(£,v), if S v > s 
lp(£v,8 — S v ) if S v < S, K v = (60) 
TO (£«)X (< vl >j s - s v) ®t v X (< u2 > ; s - Sv), else 

where the product ®j and factors m(£) are defined in lf5^|l . respectively. 

Since the underlying discrete parameter branching process defined by the 
equally likely Bernoulli offspring distribution is a critical binary Galton- Watson 
process, the recursion will terminate in a finite number of iterations with prob- 
ability one. In particular, X (9,t) is a.s. a finite product of values of \o and/or 
<p evaluated at randomly selected Fourier frequencies selected by the branching 
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<21> 



<1> 



K <i> = 
5'<i> 



S g ,K 6 = 1 



<22> 



K<2> = 1 



S<2> 



<2> 



C<i> 6 = C 



£<2> 



R 3 



Fig 1. Schematic of a tree indexed branching random walk. 



E.C. Waymire/ Probability & incompressible Navier-Stokes equations 



20 



random walk. For example, the functional evaluation of the sample tree in the 
figure is given by 

X(M) = ■m(^e)m(( <2> )(p{Ui> ,t - S <ly - S 9 ) ® ie [xo(£<2i>) ® 4<2> Xo(£<22>)] 

In general X (•) is a random multiplicative functional of scalar values m(-) and 
Fourier transformed initial data and/or forcing (vector) values over the vertices 
of a multi-type branching random walk tree tq (t) initiated from a single progen- 
itor of type £g = In general the scalar and vector valued factors are evaluated 
at the wave-number (type) of the respective vertices appearing in the tree Tg(t), 
with the initial and forcing terms appearing at the end-nodes. 

This generalizes branching random walks in the sense of LeJan and Sznitman 
jS3 for n = 3 dimensions where = |£|~ 2 ,£ G W h = R 3 \{0}. The essential 
requirement for this approach is that the following expected values exist 



The simplest proof of this uses the strong Markov property to check that if the 
indicated expected values are finite then IjFNSjl is satisfied. Existence of such ex- 
pected values has been obtained by Bhattacharya et al. |S] through constructions 
of a particular class of the Fourier multipliers, referred to as majorizing kernels, 
and defined below, from which one also obtains uniqueness in a naturally asso- 
ciated function space. Uniqueness is proven by arguments which exploit discrete 
parameter martingale structure in successive truncations by generation of the 
branching random walk product. The point is that uniqueness proofs are largely 
insensitive to the holding time distribution. In the original development of the 
theory the holding times between branchings were exponentially distributed 
with holding time parameters determined from the principal part (Laplacian) 
of the equation as above. However, recent extensions by Bhattacharya et al. [Jj, 
Orum |45j . and Ossiander |47| have been obtained based on semi-Markov cas- 
cades, which in turn have lead to new representations of both local and global 
multiplicative cascade solutions to QJ. 

A first order approach to obtain finite expected values of the branching ran- 
dom walk cascade will be seen to result from the observation that the product £g>£ 
satisfies < |u>| |z|, w, z € C", and the coefficients m(£) may be controlled 

by selecting Fourier multipliers such that m(£) < 1 or, equivalently, 



for suitable B > 0. Such h is referred to as a majorizing kernel with constant 
B. The choice of B is linked to the size of |mo| and may be standardized to 
B = 1 in this context by consideration of suitable initial data. In any case, h is 
a standardized majorizing kernel if and only if for some, and hence any, B > 
Bh is a majorizing kernel with constant B. 

The following slightly more general definition is suitable for extensions to gen- 
eralized Navier-Stokes equations with fractional Laplacian, semi-linear equations 
and for considerations of local solutions. 

Definition. A positive locally integrable function h on Wh Q R™\{0} whose 
closure Wh is a semi-group and such that 



(61) 



h*h(£)<B\Z\h{t) 1 i&W h 



(62) 
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i h is continuous on Wh 

ii h * h > a.e. on Wh 

iii h * h(£) < for £ £ Wh and some real exponent 9 and some 
B > 0, 

will be referred to as an FNS- admissible majorizing kernel with constant B and 
exponent 9. We define h = on Wfc and refer to Wh as the support of ft,. Those 
majorizing kernels which are defined and positive for all £ ^ are said to 
be fully supported. 

Formulated in these terms, the results of LeJan and Sznitman |37j may be in- 
terpreted in terms of two exponent one, standardized fully supported majorizing 
kernels, 7r 3 /l£| 2 and ae~ Q '^' /27r|£ |. These kernels are respectively non-integrable 
and integrable, with equality in (iii) of the above Definition. The former is the 
Fourier transform (in the sense of distribution) of the Riesz potential (convolu- 
tion kernel) defining the integral operator (— A) -1 , while the latter is the Fourier 
transform of the Bessel potential defining the integral operator (/ — A) -1 ; c.f. 
Folland [HO pp. 149-154]. 

Apart from their role in existence, uniqueness and expected value represen- 
tations valid in particular function spaces, the majorizing also play a role in 
tracking such structure as regularity, analyticity, support size, complcxification, 
etc. of initial data in the evolution of solutions through subspaces corresponding 
to such more detailed specifications as can be reflected in their Fourier trans- 
forms. More specifically, let us define a Banach space Th,~i,T with a norm that 
depends on a Fourier multiplier 1/h as the completion of the set 

{v G S' : t) = 0, £ G W£, \v\ T = sup < oo} (63) 

0<t<T 

under the indicated norm, where 7 G {0, 1} serves to conveniently index two dif- 
ferent norms we wish to consider. Here S' is the space of tempered distributions 
on R™. The definition of majorizing kernel and of the function spaces Th,-y,T 
imply that the product of distributions in J~h,-y,T is itself a distribution. Also, 
implicit to the definition of the Banach space Th,~i,T is the requirement that 
tempered distributions belonging to this space have Fourier transforms which 
are functions. 

In passing we wish to point out that the notion of limit space referred to 
in Section ^fia-s a natural generalization in this context. In particular if h is a 
majorizing kernel of exponent 9 = 1 then h\(£) = A -2 /i(£/A) is also a majorizing 
kernel of the same exponent. Moreover if u G Th^,T then ua G Th x , lt T and 
|u|^ t = |uA|jr h T - Thus, while an exponent one majorizing kernel h such 
that h — h\ defines a limit space !Fh,i,T m the usual sense, the norm relation 

Hw = M»w < 64 ) 

is more general. 

We will conclude this section with statements of general existence/uniqueness 
theorems possible within majorizing spaces. The main results use majorizing 
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kernels of different exponents to establish existence, uniqueness and regularity 
properties of the solutions of the l|FNS|) . Moreover these solutions have the 
expected value representation in terms of the multiplicative stochastic functional 
X {0, t) of a multitype branching random walk in Fourier wavenumber space. 
This will be followed with some basic structure theorems which are useful for 
the construction of majorizing kernels. 

The first result illustrates the use of a majorizing kernel h of exponent 1, to 
obtain existence and uniqueness of the solution given by the stochastic cascade 
representation for small enough initial data and forcing on a time interval that 
is solely constrained by the length of time for which the forcing remains small. 
Specifically one has the following theorem. 

In the statements of these results, (— A) _/3 , < (3 < 3/2 denotes the negative 
power of the Laplacian defined as the integral operator with symbol |£|~ 2/3 ; i.e. 
the convolution with a Riesz potential. 

Theorem 3.1. Let h(£) be a standard majorizing kernel with exponent 9 = 1. 
Fix < T < +oo. Suppose that \u \r h T < (V2n) 3 v/2 and K-A)" 1 ^ g t < 

(\/27r) 3 // 2 /4. Then there is a unique solution u in the ball Bq(0, R) centered at 
of radius R = (v / 27r) 3 ^/2 in the space Fh.o.T- Moreover the Fourier transform 
of the solution is given by = h(^)E^ e= ^){ € Wh- 

As already noted the existence question is made trivial by the stochastic iter- 
ation so long as the indicated expected values are defined. A proof of uniqueness 
can be based on a discrete parameter martingale construction by truncated ge- 
neology adapted from LeJan and Sznitman |37j which is used to obtain a bound 
on the difference between two solutions in this space by 2P(T > n) for all 
n = 0, 1, 2, . . ., where T denotes the time to extinction of the underlying dis- 
crete parameter critical binary branching process. 

It should be remarked that some regularity properties of the solutions can 
be inferred from the particular majorizing kernel being used via Payley- Wiener 
theory, Sobolev embedding, etc. For example, note that the majorizing kernel 
^o(C) = 7r3 /IC| 2 gives existence and uniqueness, but no control over tracking reg- 
ularity properties of the solutions from those of initial data. However, solutions 
obtained using the majorizing kernels h^ = \^f- 2 e~ a ^\o < /3 < l,a > 0, 
maintain the C°°— regularity of the initial data, as can be seen from the bound 
on the Fourier transform of the solution Moreover one may construct, for exam- 
ple, smooth compactly supported initial data whose evolution is more efficientlly 
dominated by a kernel h^ with (3 < 1, whose decay at infinity is consistent with 
the compactness of support property; see Bhattacharya and Rao pp. 83-89]. 

On the other hand, as will be elaborated below, working in the function 
spaces Th.x.T it is possible to use majorizing kernels to obtain spatial analytic- 
ity of the solution. However, it should be noted that the size constraints imposed 
on the initial data and forcing to obtain analyticity are substantially more re- 
strictive than those required in Theorem 13.11 i.e. existence and uniqueness of 
such solutions is within a smaller ball of the function space. Specifically one has 
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Theorem 3.2. Let h(£) be a standard majorizing kernel with exponent 9 — 
1. Fix < T < +oo. Assume {e^uoix)^ 1 T < (v/ f )3 pW 1/2l/ and that 
\(-A)- 1 g(x,t)\ :Fh < { ^ )3 pv 2 e~ 1/2u for some < p < 1. Then there is a 
unique solutionu in the ball Z?i(0, R) centered atO ofradiusR= (p/2)(\/27r) 3 ^x 
e _ 37; in the space Fh,i,T- 

Under the conditions of Theorem 13.21 the asserted solution satisfies the fol- 
lowing decay condition 

e^\H(£,t)\ 

sup sup — — < oo (65) 

0<t<T£6R3 /i(?) 

Thus Theorem 13.21 provides another approach generalizing that of Lemarie- 
Rieusset [35], Foias and Temam [3D] to obtain conditions for regularity in the 
stronger form of spatial analyticity. In particular, by means of the following key 
inequality along the lines of Foias and Temam |2U] and Lemarie-Rieusset |35) . 
for < s < t,£,T) E R 3 , 

CX P {-!AS|£| 2 - VT~S\£ - 7/| - Vt~ S \ V \} < eW e -V*|el e -^|er/2 (gg) 

the majorization of u(£, t) can be shifted to e v/ *'^'u(^, t); see Bhattacharya et al. 
6, Proposition 4.1, p. 5031]. As a consequence, for example, if e - *'^ft(£) e L 1 
for some 6 then, for times \ft > S, u(x + iy, t) is complex analytic in the region 
x + iyeC 3 such that |y| < — 5. Specifically, 

u(x + iy,t) = / e J ( x+jy >«u(£,i)d£ 

^e-^W+^e^uit, t)dt (67) 



'R 3 

is a complex analytic extension since 

< e- s ^h(£) (68) 



for |y| < — S. In particular, in the case e _l5 ^'ft.(^) 6 L 1 for all 5 > 0, 
analyticity holds for all times t > 0. This is the case for h(£) — 1/|£| 2 treated by 
Lemarie-Rieusset jSSJ) as well as for such majorizing kernels larger at infinity as 
indicated in examples below; c.f. Example (viii). 

Two additional corollaries to the theory for global solutions follow imme- 
diately from the stochastic cascade representation: (a) time-asymptotic solu- 
tions, c.f. Bhattacharya et al. |7j, Cannone [S] and references therein, and (b) 
self-similar solutions, c.f. LeJan and Sznitman [3Tj . Cannone |H] and references 
therein. Specifically, for the case of (a), under the majorizing conditions of Theo- 
rem ^. II with T = oo, suppose further that lim^oo g(^ t) = <?oo(£) exists for each 
£ =/= 0. Then, since the underlying discrete parameter binary branching is criti- 
cal, lim^oo X (0, t) exists finite random product. Moreover, under the 
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majorizing conditions one has for each t > 0, with probability one |X {0,t)\ < 1. 
Thus, by Lebesgue's Dominated Convergence Theorem 

Xoo(0 := lim x(€,t) 

t — >oo 

exists for each £. Now again apply Dominated Convergence to (FNS)h to obtain 

}ds. 
(69) 



Multiplication through by h(£) establishes that 

t — ►OQ 

exists and satisfies the steady state Navier-Stokes 

fioo(f)=/ e-^l 2s {|e|(27r)-i / Uoo(v)^oo(A-v)dv+9oo(0}ds. (FNS)^ 
Jo Jn 3 

With regard to (b), mild self-similar solutions to Navier-Stokes of the form 

, , 1 , . _L - * 

u{x,t) = Vt { Vt P(X '^I P( 7! ) 

were originally considered by Leray in connection with backward-in-time 
evolution of Navier-Stokes as an unrealized possible approach to determine sin- 
gular flow. Cannone :9; provides a brief historical overview and ultimate demise 
of backward self-similarity, but the existence of forward self-similarity. In this 
regard the simplicity of the branching random walk cascade is overwhelming. 
Specifically, one observes that for majorizing kernels h defined by equality in 
<|{j2l) satisfying the limit-space scale invariance 

MO = A- 2 M£/A), A > 0,^0, (70) 

it follows that (i.) m(£) is constant and therefore trivially invariant under scale 
change £ — > £/A, A > 0; (ii)a®£ b depends on £ only through eg = t|t ; under the 
scale change — > £„/A, A > 0, the distribution of the holding time S v changes to 
that of X 2 S V ; and, finally, the spatial transition probability ^j^jj^y^dri is invari- 
ant under re-scaling of Fourier frequencies. Thus taking initial data 
and forcing g(£)/h(£) (constant in time) to both be homogeneous of degree 
one obtains that 

u(£,i) =A- 2 u(A- 1 £,A 2 t). (71) 

Cannone [2] notes the heuristic connection between time-asymptotics and self- 
similar solutions obtained by observing that since for a global solution u,p, 
the space-time rescaling u^ A '(x, t) := Au(Ax, X 2 t), p( A )|(x) = A 2 p(Ax), is also a 
solution, it follows that u(°°)(x, t) := limA->oo u( A )(x, t) is self-similar whenever 
the indicated limit exists. In particular limt_, (X > 1 3 u(t%x, t) = u(°°)(x, 1). 

Finally let us record local existence and uniqueness from more relaxed condi- 
tions on the majorizing kernels as illustrated by the following result; c.f. Bhat- 
tacharya et al. [E]. 
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Theorem 3.3. Let h(£) be a standard majorizing kernel with exponent 9 < 1. 
Fix < T < +00,7 G {0,1}. Assume e vtA Uo(x) G Th^,T & n d for some 1 < 
(3 < 2, (— A) _ § g[x, t) G J~h,~f,T Then there is a < T* < T /or which one has 
a unique solution u G J-h,j,T, ■ 

Another variation on the general approach presented here involving semi- 
Markov branching, in the sense that the holding time distributions are not re- 
quired to be exponentially distributed, leads to conditions for a local existence 
and uniqueness theory in all dimensions.; see Bhattacharya et al. 0, Orum |45|. 
Ossiander g7J. 

Some sense of the class of majorizing kernels may be derived by noting from 
Holder's inequality that the set of fully supported majorizing kernels with a 
given exponent is a log-convex set. Also if is a majorizing kernel then so 
is ce a '^n(£) for arbitrary fixed vector a and positive scalar c. One may also 
show that the only fully supported homogeneous majorizing kernels in n > 3 
dimensions are those of degree n — 1. Along these lines the Theorem 13 .41 below 
may be viewed as something of a "toolkit" for constructing examples. 

The family of standard majorizing kernels of exponent 6 on R" will be de- 
noted by 

H n , e = {h : W h - (0, oo) : h * MO < for all £ e W h C R™}. 

Theorem 3.4. [Majorizing Kernel Toolkit] 

1. Suppose that {qj : 1 < j < m} satisfies qj > 0, and Qj = 1- F° r 
hj G Hnfij, j = 1, -,m, 

h(0= Rhf® 

with support Wh = nJLiWh r 

2. For hj G Tinfi: j = 1, —, m, 

M0= n hf(0 eH n .e 

wii/i support Wh — (^"LiWhj ■ 

3. Fix n > 1. Suppose that ki , . . . fc TO is a partition of n and for each j = 
1, . . . , m, /ij is in Tiktfij ■ Then 

m j 

HO = II g = (6, • • . U) for e R 

J=l 

is in Hn^f^Bj with W h = W hl x ... x 
^. //A is a n x n invertible matrix and h G W n ,0, i/ien defining \\A\\ = 
sup{|Ax| : |x| = 1}, 

h A (0 := \detA\ ■ \\A\\- e h(M) G W„,e 
luii/i support Wy lA = {A^ 1 ^ : £ G W/J. 
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5. If h G H n .o and ip : R™ — > [0, oo) satisfies ip(£) < '0( 7 7) + — 7 ?) / or a ^ 
77, £ € W/,, tften 

6. J/fte H n ,e, fen e Q «/i(0 G H nfi for any fixed a G R" , 

7. If h G 7i n ,e i/ien /or any pseudo-metric p on a subset of R 3 contain- 
ing Wh, e~ ap ^°'£)/i(£) G H n ,e /or any a > and £o fixed. In particular 
e - a \i\ h(O £ Hnfi for any a> and < (3 < 1 

In general, existence and classification of majorizing kernels is non-trivial. 
For example, it can be shown that any piecewise continuous h G Hi t i must 
have Wh = (0, oo) or Wh = (— oo,0). This illustrates the tradeoff between 
n and 8; if exponent 8 > 0, the existence of majorizing kernels with support 
R"\{0} is problematic for smaller values of n. There are however fully supported 
majorizing kernels of exponent 8 — for all n > 1. Using the toolkit Theorem ^. 41 
one may identify a number of majorizing kernels; see Bhattacharya et al. joj. 

Examples [Majorizing Kernels]. 

1 MO = ihrfr+p) e n ho with W hl = R. 

ii For n > 1, ftn(0 = (2tt)-" n" =i (l + £f )~ 1 G ft„,o with W hn = R". 
in MO = „+i r( ^ } ,+i G H n , 

iv For £ ^ 0, < /? < 1, a > 0, ftj^O = \£f- 3 e- a W G W 3 ,i with 
W = R 3 \{0}- 

v For each 8 G (0, 1), < /3 < 1, and a > 0,£ 7^ 0, one has that hi a 2(£) = 



|^|8((3-2) „-c,8|5|' 



(2^)3(i-»)n, a = i(i+c|) (i - e) 



G H 3 ,e with FF = R 3 \{0}. 



vi For each 8 G (0, 1), < (3 < 1, and a > 0,£ 7^ 0, one has that h£*2(£) 



(i+iei 2 ) 2(1 - 



G K 3 ,e with VF = R 3 \{0}. 
vii For n > 3 and 7) with < (3 < 1 and 1 < 7 < 1 + (3, one has that 



V/3,7(0 = /t > o* Vl ~ 1 e~ t ''~ ICI " /t * G W„,i with Wfc = R"\{0} 
viii Let G be defined a.e. on the unit sphere with G(8) — > 00 as approaches 
the points (0, 0, ±1), (0, ±1, 0), (±1, 0, 0), respectively. Then one has that 

The families of non-fully supported non-radial kernels illustrated in Example 
(viii) are larger kernels that permit broader existence and uniqueness results for 
given initial data u of (|FNS(I . In particular the growth of hs,i along particular 
directions is much larger than ft.o(0 = VIO 2 - Transforming hgi via a rotation 
permits such growth in any direction. 

The kernels in Example (vii) are closely related to the Bessel kernels of Aron- 
szajn and Smith 2 . They can also be combined with other kernels of the Ex- 
ample (vii) along the lines of the previous proposition to construct kernels in 
Hn,e for < 8 < 1. One may apply the Laplace method for estimating integrals 
to show that the Bessel type kernels h = h n ^ tl are also regularizing kernels in 
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the sense that the distributions in the corresponding function space !Fhfi,T are 
C°°— functions. 

In view of the role of majorizing kernels in providing bounds on the Fourier 
transformed forcing and/or initial data, the theory contains a dual problem 
which is to identify classes of divergence free vector fields in physical space 
which are so dominated. This section is concluded with some examples in this 
category. 

The first example is a class of solenoidal vector fields on R 3 whose Fourier 
transforms are dominated by /i3,/3, 7 (0- 

Example. Fix < (3 < 1 and 1 < 7 < 1 + f3. For 1 < j < 3 let mj(t) be measur- 
able functions on [0, 00) such that |mj(i)| < ti~ 1 e~ tf and J t>0 t~ 3/>2 \mj(t)\dt < 
00. Let v(x) be the vector field whose components Vj(x) are the Laplace trans- 
forms of rrij(t) evaluated at |x| 2 /4; that is 

/•OO 

Vj (x) = / e- tlxl2/4 mj (t)dt 
Jo 

Let u(x) be the divergence free projection of v(x). Then 

|«(OI < Cfc3,/J, 7 (0- 

For the next example we consider majorization by the kernels hfj^ . 
Example. Let M. denote the space of finite signed measures on R 3 with total 
variation norm || ||. Let < (3 < 1 and denote the "Fourier transformed Bessel 
kernel" of order (3 by Gf3(x) = (1 + |a;| 2 ) _ ^. Then for each g = Gp * [i, fi e M, 
one has for (3 = 1, a e (0, 1) and for (3 e (0, 1), a > 0, 

for a constant > 0. In particular, if v <E L 1 is a divergence free vector field 
then g = Gp*v is also a divergence free vector field whose Fourier transform is 
dominated by h^"\ 

Finally the following example uses the majorizing kernels to provide 
smooth solenoidal vector fields, including some with compact support, for which 
the theory applies. 

Example. Let irij(t),t > 0,j = 1,2,3, be measurable functions such that 
J °° e~\ x \ 2t \m,j(t)\dt < oo,x e R 3 ,j = 1,2,3. Define a vector field with com- 
ponents Vj,j — 1,2,3, by 

/•OO 

Vj(x) = / e^ xl2t m: j(t)dt,x e R 3 . 
Jo 

Let u be the divergence free projection of v. Then, 

(i.) If \mj{t)\ < ct-i then \v,j(£)\ < c'/i a) (0 for some c' > 0,j = 1,2,3. 
(ii.) If \mj(t)\ < ce- 2 " 2 * then |%(£)| < c'h ( *\£) for some d > 0, j = 1,2,3. 
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(iii.) For arbitrary e > there is a smooth probability density function k e 
supported on [— e,e] 3 such that 

\k(0\ <c(/3,e)exp{-|eC| /3 },eeR 3 ,c(/3,e) > 0. 

Let v be any divergence-free integrable vector field such that |w(£)| < c|£| -2 , £ ^ 
0. Then the component-wise perturbation u — k e *v is a divergence-free infinitely 
differentiable vector field such that < c'h^l^), for a = e 13 and some 

c'>0,i = 1,2,3. 

4. Some Concluding Remarks 

We wish to conclude with pointers to two major problem areas pertaining to 
incompressible Navier-Stokes equations which also fall within the subject area 
of probability theory. 

Probabilistic methods in the study of Navier-Stokes equations are not new, 
howeverr new ideas have emerged recently which are based on objects of broad 
probabilistic interest such as multi-type branching random walks. The theory 
to date is in a form closely tied to the contraction map methods of Kato |30|. 
so much so that most can be obtained by either route. Bounding the stochastic 
cascade product by one is a bit simpler than the contraction argument only in 
the case that the contraction is not strict. Similarly, as illustrated in this survey, 
owing to the a.s. finiteness of critical branching, uniqueness, time-asymptotic 
limits and certain self-similarity structure are extremely simple observations 
within the stochastic cascade framework. 

The iterates of the Picard scheme can be probabilistically identified as ex- 
pectations over truncations ( by generation) of the branching random walk. In 
particular, the n— th iterate is represented by an expected value restricted to 
averages over trees whose evolution by time t consists of only vertices v such 
that \v\ < n and such that the n-th generation vertices |u| = n survive beyond 
time t, i.e. Y^j=a &v\j > see Proposition 4.3 of Bhattacharya et al. jf) . This 
and the uniqueness proof are two instances in which minimal but significant 
consideration is given to this underlying branching evolution. 

However suitable exploitation of the genealogy of the branching random walk 
in ways which take fuller advantage of incompressibility has the potential to 
yield significant new results. Problems of this latter type are indeed at the next 
frontier of this research. 

The Fourier representations are naturally constrained to free space and/or 
periodic boundary conditions. A second major problem for probabilists is the 
continued quest for representations in physical space which are amenable to the 
imposition of boundary conditions of various types. 
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